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equilibrium, or nonequilibrium. Also the effects of the diffusion 
model, surface catalyticity , and mass injection on surface transport 
and flow parameters are considered. The equilibrium calculations for 
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brium calculations where the diffusion is multicomponent. 

The viscous shock-layer equations are solved by using an implicit 
finite-difference scheme. The flow is treated as a mixture of inert 
and thermally perfect species. Also the flow is assumed to be in 
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All calculations are for a 45° hyperboloid. The flight conditions 
are those for various altitudes and velocities in the earth's atmos- 
phere. Data are presented showing the effects of the chemical models, 
diffusion models; surface catalyticity; and mass injection of air, 
water, and ablation products on heat transfer; skin friction; shock 
stand-off distance; wall pressure distribution; and tangential velocity, 
temperature, and species profiles. 
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I . INTRODUCTION 


Analysis of the flow about a hypersonic vehicle must account for 
the Interactions that occur between the reacting outer flow and the 
vehicle surface. This problem is complicated in many cases because the 
vehicle surface is either an ablative or transpiration-cooled surface. 

In such instances, gaseous species that may have a wide spectrum of 
molecular weights are injected into the external flow. Consequently, 
assumptions concerning the chemistry and diffusion models may signifi- 
cantly affect the accuracy of an analysis for such a problem. 

Numerical solutions to the aforementioned problem have been either 
that of the inviscid-boundary-layer solutions (refs. 1 through 6, for 
example) or the viscous shock- layer solution (refs. 6 through 16, far 
example). The viscous shock-layer solution as originally suggested 
by Cheng (ref. 16) and later modified by Davis (ref. 17) Is a most 
appealing method because one set of equations uniformly valid through- 
out the shock layer is used. Furthermore, the viscous-inviscld inter- 
actions are accounted for in a straightforward manner with or without 
mass injectloi-. 

The most recent numerical solutions of the viscous shock- layer 
equations have been those of Davis (ref. 8), Whitehead (ref. 9), 
Dellinger (ref. 10), Blottner (ref. 11), and Adams (ref. 6). References 
8 and 9 used the most extensive shock- layer equations since the solu- 
tions are not restricted to the thin shock-layer approximations as in 
the other stated references. The solutions of references 8, 9, and 10 


1 
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are for the stagnation and downstream flow, whereas those of references 
6 and 11 are for the stagnation streamline only. However, the above 
viscous shock-layer solutions considered rather simple chemical systems 
of species. For example, reference 8 considered atomic and molecular 
oxygen species using finite-rate chemistry; reference 9 considered the 
effects of injecting argon, air, and helium into air where air was 
treated as one species and the chemistry was either frozen or equili- 
brium; reference 10 considered nonequilibrium air as a reacting mixture 
of seven chemical species with constant but arbitrary Prandtl and Lewis 
numbers; and references 6 and 11, which presented stagnation solutions, 
considered the same chemical species as reference 10, but reference 
6, in addition, considered multicomponent diffusion and the effects of 
inert gas injection. Yet, none of the above viscous shock-layer 
solutions for either equilibrium or nonequilibrium chemistry has all 
of the following features: downstream solution capability, mass 

injection, multicomponent diffusion, and second-order terms in the 
normal momentum equation. However, this study presents numerical 
solutions of the viscous shock-layer equations where all the above 
features are accounted for. The chemistry is treated as being either 
frozen, equilibrium, or nonequilibrium. The equilibrium air calcula- 
tions that treat the diffusion as multicomponent provide solutions 
previously unavailable. The same can also be said for the downstream 
nonequilibrium air calculations where the diffusion is multicomponent. 

The viscous shock-layer equations and the implicit finite-difference 
scheme used to solve the equations are essentially those given by 
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Davis in reference 8. The flow is treated as a mixture of inert and 
thermally perfect species. Also, the flow is assumed to be in vibra- 
tional equilibrium. 

All calculations are for a 45“ hyperboloid. The flight conditions 
are those for various altitudes and velocities in the earth's atmos- 
phere. Data are presented showing the effects of the chemical models; 
diffusion models; surface catalyticity ; and mass injection of air, 
water, and ablation products on heat transfer; skin friction; shock 
stand-off distance; wall pressure distribution; and .tangential velocity, 
temperature, and species profiles. 
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II. DEVELOPMENT OF THE GOVERNING EQUATIONS 

The conservation equations and their associated boundary condi- 
tions are presented for the laminar viscous shock layer about an hxis- 
ymmetric or two-dimensional body. These equations account for mass 
injection and treat the diffusion as either multicomponent or binary. 
Equations convenient for solving frozen, equilibrium, and nonequilibrium 
flow are presented. Also, the non-dimensionalizing quantities are 
given, along with the transformation used to facilitate the numerical 
solution of the equations. 

Flow Field Equations 

The conservation equations that describe a reacting multicomponent 
gas mixture can be found, for example, in reference 18 or 19. The 
viscous shock-layer equations as used herein are obtained by expressing 
the conservation equations in the body— oriented coordinate system shown 
in figure 1. Then these equations are nondimensionalized in each of 
two flow regions. That is, the equations are nondimensionalized with 
variables which are of order one in the region near the body surface 
(boundary layer for large Reynolds numbers), and then the same set of 
equations is nondimensionalized with variables which are of order one 
in the essentially Inviscid region outside the boundary layer. Terms 
in each of the two resulting sets of equations are retained up to 
second order in the inverse square root of a Reynolds number. By 
combining these tx^o sets of equations such that terms up to second 
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order in both the inner and outer regions are retained, a set of equa- 
tions uniformly valid to second order in the entire shock layer is 
obtained. 

The resulting equations are of a hyperbolic-parabolic nature. 
However, if the thin-shock-layer approximation is made (an approximation 
to the normal momentum equation), the equations are parabolic. The 
parabolic equations can be solved by using numerical methods similar to 
those used in solving boundary- layer problems. Then, an iteration 
procedure is used to remove the thln-shock-layer approximation, as will 
be discussed in a later section. 

The derivation described above and the subsequent equations are 
the same as given by Davis (ref. 8) except that the energy equation, 
in addition to being formulated in terms of temperature, is formulated 
in terms of total enthalpy and the species equations are also formu- 
lated as elemental equations. The viscous shock- layer equations for a 
chemically reacting multicomponent mixture are: 

Global Continuity 

[(r + ncos0)^pu] + [(1 + nic)(r + ncos6)^pv] = 0, (2.1) 

os on 

s-momentum 

[ u 3u + uvK ~ l , 1 3P _ 

1+nK 9s l+nKj- 1+nK 9s 
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n-momentum 
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which becomes 


(2.3b) 

3 n 1 +nK 


if the thin-shock-layer approximation is made, 


Energy (temperature) 
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Species Continuity 
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h ({^ 


Where is the diffusion mass flux of species i, and 

State 



(2.5) 



The above equations are in a form that is convenient to apply when 
the flow-field chemistry is nonequilibrium. In the energy equation and 
the species continuity equations, the term appears. This term 

represents the rate of production of species i due to chemical 
reactions. For frozen flow the rate of production terms is zero. 
However, when the gas is assumed to be in chemical equilibrium, the 
production terms cannot be obtained from the chemical kinetics. Thus, 
the production terms in the energy equation are eliminated by formulat- 
ing this equation in terms of enthalpy. 

The energy equation in terms of total enthalpy becomes 


r u 8H 3 h"| 

1 1+nK 3s 3n 1 
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where the total enthalpy as used in the above equation is given by 



( 2 . 8 ) 


This particular formulation of the energy equation reduces the coupling 
between the energy and species continuity equations (ref. 20) over 
that of equation (2A) . 

The production terms that appear in the species continuity 
equations are eliminated by Introducing the concept of elemental mass 
fractions as expressed by Lees in reference 21. As long as no nuclear 
reactions occur, the elemental mass fractions remain fixed and unchanged 
during chemical reactions. The relation between the elemental and 
species mass fractions is given by 


NS 




^5, ^i 


(2.9) 


where 

<S . „ is the number of atoms of the Jlth element in species i, 

Xa/ 

is the molecular weight of the 2-th element, 
is the molecular weight of the ith species , 
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and 

is the mass fraction of species i. 

The elemental continuity equations for the elements can be obtained by 
multiplying equation (25) by summing over i. The 

resulting elemental continuity equations are. 

P|_l4tiK 8s -"^an J = (i+nK)(r+ncos0)J 

^{l + nK} (r + ncos6pJj^)j , (2.10) 

where 



NS M* 

'^i • 


( 2 . 11 ) 


Introducing the elemental mass fraction not only eliminates the pro- 
duction terms in equations (2.4) and (2.5) but also reduces the number 
of equations to be solved, since there is one equation (eq. (2.10)) 
for each element rather than one equation (eq. (2.5)) for each species. 

In addition to the above system of equations, additional equations 
are required to specify the mixture thermodynamic quantities and the 
diffusion mass flux quantities. The thermodynamic relations are as 


follows : 
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Mass Fractions 


C. = p^/p 


Frozen Specific Heat 


C = 
P 


NS 


i=l 


C.C , 
i pi 




Molecular Weight 


NS C. 


and 


Enthalpy 


NS 

h = E h.C^ , 
i=l " ^ 


( 2 . 12 ) 


( 2 . 13 ) 


( 2 . 14 ) 


( 2 . 15 ) 


where h^ includes the enthalpy of formation of the ith species. 

The mass flux relative to the mass average velocity when thermal and 
pressure diffusion are neglected can be written in terms of the multi- 
component Lewis numbers, j » 22) as 


NS M*_,M*. 9x. 


j = l M*^ 


¥± 


( 2 . 16 ) 
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or as given in reference 1 as 


NS 3C, 

J. = - ^ E Ab . 

1 Pr , T ik 3n 
k=l 


(2.17a) 


where 


Ab 


ik 


[ M*. 

M* 


M*. 

+ U - M. 


A os -j 

i E L. .C. 
k/j = l U 


for i = k 


for i 7 ^ k , 


(2.17b) 


and 


NS 

E C./M*. 

j=l 


2 2 




NS 

E C./C..M*, 
1 J iJ j 


(2.17c) 


and C.. are the binary Lewis numbers. The relative mass flux for 
ij 

the elements can be written as 


h- 


f B 

Pr L J ’ 


(2.18a) 


where 


NS M*. ^ 

\k " ^ iS, ^^ik ’ 


(2.18b) 
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L. - L for i = k 

1 

FM*. / M*.\ NS ”1 

and L is an arbitrary constant. When the diffusional mass flux for 
the elements is written as indicated, the elemental continuity 
equation can be arranged in a convenient form for numerical solution. 
Equations (2.17) and (2.18) can be simplified to 




(2.19) 


and 




( 2 . 20 ) 


respectively, whai the binary Lewis numbers, > are all equal to 

^ or an "effective" Lewis number is used that is equal to 

The previous equations have been nondimens ionalized by using the 
following relations ; 


u* = uU* 

OO 


> 


(2.21a) 


V* = vU* , 

OO ' 


(2.21b) 
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T* = TU* /C* , 

“ P 

oo 


P* = pp *u* , 

^ ^ ^OQ OO * 


P* = PPco* 


y* 


yy* 


fu* ^/c* 

l OO p 




K* = Ky* ^C* 
ref p 

CO 


J 


c* = c c* , 

p p 


h* = hU* 


= \ P*J^*J^* > 


J*. = j.y* ./a*, 
i 1 ref 


s* = sa* , 



K* = Ka* , 


(2.21c) 

(2.21d) 

(2.2ie) 

(2.21f) 

(2.21g) 

(2.21h) 
(2.211) 
(2.21j) 
(2.21k) 
(2.211) 
(2.21m) 
(2. 2 In) 


O 
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and 

r* = ra* . (2.21o) 

Also, four dimensionless parameters appear in the shock- layer equations. 
They are 


and 


Pr = C* y*/K* , 

P 


e 



ref 


0 * U* a* 

^ CO CO 



P*C*^*, ,/K* 
P ij 


» 


(2.22a) 


(2.22b) 


= p*C *D* ./K* (2.22c) 

ij P ij 


To simplify the numerical computations, a transformation is applied 
to the foregoing equations. This is accomplished by normalizing the 
independent and dependent variables with their local shock values. 

When the normal coordinate is normalized with respect to the local 
shock stand-off distance, a constant number of finite-difference grid 
points between the body and shock are used. Also, the need for inter- 
polating to determine shock shape and the addition of grid points in 
the normal direction is eliminated. 

The transformed independent and dependent variables are 
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n = n/n^ 


C = s , 


U = u/u^ , 

V = v/Vg , 

p = P/Pg . 

p = p/p, , 

s 

T = T/T , 
s 

H = H/H 

s 

ij = y/Pg , 

K = K/Kg , 


(2.23a) 

(2.23b) 

(2.23c) 

(2.23d) 

(2.23e) 

(2.23f) 

(2.23g) 

(2.23h) 

(2.231) 

(2.23j) 


and 


C = C /C 
P P P. 


(2.23k) 


The transformations relating the differential quantities are 
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3s ~ Hg 3ri 


^ ^ 1 ^ 

3n 9ri ’ 


and 
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2 2 
3n" 


where 

dn 

I s 

n = -Jc~ 
s d^ 


(2.24a) 


(2.24b) 


(2.24c) 


(2.24d) 


The s-momentum, energy (temperature and enthalpy), species continuity, 
and elemental continuity equations (eqs. (2.2), (2.4), (2.7), (2.5) and 
(2. 10), when written in the transformed 0 plane, can be put into 
the following standard form for a parabolic partial differential 
equation: 


2 

3 W 

2 

9n 


“i 3n “2'^ 


^ _ _ 
“3 "^“43^ ® 


(2.25) 


where W equals u in the s-momentum equation, f in the temperature 
energy equation, H in the enthalpy energy equation, in the 

species continuity equations, and in the elemental continuity 

equation. The coefficients through may be written as follows: 
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s -momentum 


1 n K 

1 3y , s 

a, = - + 


+ 


in COS0 
s 


1 - 9n 1+n r|K r+n ricos6 

y s s 


n D u n' — n p V — 

+ s s s s pun s s s pv 

e^yg(l + n^nx) y ij 


(2.26a) 


Kn 


“2 ■ ■ ( 1 +ngXn)- 3 n ■ ) 2 " (ri^gncosO) (H^gHK) 

N \ S ' 


s 1 3y 


2 2 
K n 


jcos9n K 
s 


2 , 
p n u' 
s s s 




iip 

y 


n p V X 
s s s 


2 

e y 


£Z 

y 


^3 " ■ 


P n^ 
s s 


■=\“s " 


Fsp 3 P ~1 

L9? Pg ■ 9nJ ’ 


“4 ■ _2 


s 




PH 

y 


(2.26b) 


(2.26c) 


(2.26d) 


Energy (temperature) 


NS 


a. 


1 9K r 

K 

. jcos0 

n ,E, J.C 
s i=l 1 p^ 


= ^ 3n ^ "s L 

. 1+n^nx 

r+n T 1 COS 0 

s • 

K K 
s 


n p C pC 
s s p p 

s 

1 V V - 

L s 

V's^l 


(2.27a) 

2 - 
e K K 
s 

1+n r|K J ’ 



f 

a.T' 

_ 4 s 

2 . 

^ s^s2 



(2.27b) 


e K K 
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a, = - 


2 . T 

n w T W-, 
s si 1 , 

2 - , 
n VV 
s s 


2 


L ^ 

e T K K 

T K K 

L- s 

s s 

s s 


n^ u P u 

r 

■p’ ' 

s s s 


s 


K K (1+n riK) 
s s s 


3P 

9 ^ 


n 9ri 
s 


] 


n V P ^5 
s s s - oP 

ATi''"’! 

s s 


(?,.27c) 


“4 


n p C u 
s^s Ps s 

e^K (l+n n< 
s \ s 


pC u 

P ■ 

K 


(2.27d) 


Energy (enthalpy) 


1 3y , 1 9Pr , 

“i = -9n-if9n 




jcos9 


l+n n< r+n Ticos0 
s s 


] 


p Pr 
s s 


“2 H 


2 

Vs 


Pr n Fn' u rjup 

s I s s 

L l+n r|K 

II ^ c 


- V pv 
s 


il 


(2.28a) 


(2.28b) 


as = 


Pr Pr 
s 




r 1 

S 

n 

_ s 


9y 


.i cos 6 


l+n r|K r+n pcosB, 
s s y 


(2.28c) 


“4 = - 


Pr n p u Prup 
s s s s 

e^y ( l+n pk) y 


(2.28d) 


where 
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¥ = 


n Pr 
s s 


[ 


NS 9C. 

E h. ^ + 


Pr i=l 


‘i 9n 


2 — 
u yu 


Pr 


Pr Pr 
V s 


] 


NS 

- E h.J. 
i=l " ^ 


2 —2 
y u Kyu 
s s 

l+n m 
s 


(2.28e) 


The above relations are for the thin shock-layer approximation. When 
equation (2.3a) is used for the n-momentum equation, the following terms 


must be added to OL^' 


Pr V n p P 
s s s s , 

2 

e y H 
s s 


rpv p u^u 

l+n r 
u l— s 


+ \ (2.28f) 

„ r|K \ s sdq n oT] I n on. 
yi_s \ s / s 


:] 


Species continuity 


1 9PL ^ 
“l “ PL 9n 


s 1 


K , J_cos 9 


l+n PK r+n ncos0 
s s 


9 


p V n pv n pun' pun 
sss^ sss s 


= - 


e^PL 

2 .1 
n w . 
s 1 


e'^PL ( l+n nK ) 


2 

e PL 


ao = 


1 r 9PM . ..../ K , — . 

PL I 9n ( l+n nK r+n ncos0/ 2 

I— ' s s 


.0 2 

jcos0 \ ^ i 




(2.29a) 


(2.29b) 


(2.29 c) 


a, = - 


n p u p .u 
s s s 

e^PL (l+n n<) 


(2.29d) 


where 
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PL = 


Pr Pr 
s 


and 


PM = 


V 


NS 
Z Ab. 


9C, 


Pr Pr k=l 


i,k 9n 


for multicomponent diffusion and 


PL = 

Pr Pr 
s 


and 


PM = 0 


(2.29e) 


(2.29f) 


(2.29g) 


(2.29h) 


for binary diffusion. 
Elemental Continuity 


= r K ^ ,jcos6 

°^1 9n 1+n r\K r+n ncosO 

PL L s s _J 


pvnpv npun’ pun 
*^s s s^ , s'^s s s 


2^^ 2'^ I \ 

e PL e PL 1 + n^nKj 


a2 = 0 


a. 




9n 


PM + 


n PM 
s 


( 


1+n nic 
s 


+ 


_i cos 6 


r+n r)cos0 
s 


)] 


(2.30a) 


(2.30b) 


(2.30c) 



21 


P u -- 

s s s £u 

“ .2 , , . \ PL 


e ( 1+n n< ) 


where 


PL = 


Pr Pr 
s 


(2.30d) 


(2.30e) 


and 


^ NS M* NS 9C 

1=1 i k=l 




for multicomponent diffusion and 


(2.30f) 


and 


PL = 

Pr Pr 
s 


PM = 0 


(2.30g) 


(2.30h) 


for binary diffusion. 

The remaining equations may be written as follows; 


Global Continuity 


9C 



r + 


n^pcosS^ 


j 


p u pu 
s s 



r + n ncos9 
s 




V pv - n' np u pu 
s s s s 



0 , 


(2.31) 
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n-momentum 


pu s - , 8v 

1+n r|ic V oC 

s _ s 


tl T| rv \ 

s ov ^ s _p_v ^ 

n u n 8ri 

s U s s 



—2 

pu 


P 

s 


p u n V 
s s s s 



(2.32a) 


which becomes 


9n 


n p u K 
s s s 


—2 

pu 


(2.32b) 


if the thin shock-layer approximation is made, and 


State 



M* 

s 


(2.33) 


The terms and w^ which appear in the temperature form of 

the energy equation and the terms ^ ^ which appear in the 

species continuity equations are quantities that involve the rate of 
production terms w^. As discussed in reference 1, 8, and 11, the 
way the production terms are written is very important to achieve 
convergence of the iteration procedure. Consequently, for the energy 
equation, the production terms were written as given in reference 8 


as 
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where k denotes the iteration for which the solution is known and 

k + 1 the iteration for which a solution is required. Accordingly, the 
NS 

term E h.w. which appears in equation (2.4) was written as follows: 
i=l 


NS 

E 

i=l 


h.w, 
1 1 


'"sl'^1 "■ 


Vs2'^2 


(2.35) 


and appears in equations (2.27b) and (2.27c). As for the species 
continuity equations, the production term was written as 



- .1 
C.w , 
1 i 


(2.36) 


and appears in equations (2.29b) and (2.29c). Hence, equations (2.34) 
and (2.36) express the production terms in terms of temperature for 
the energy equation and species mass fraction for the species equations. 
Blottner in reference 1 comments on the necessity of using expressions 
like equations (2.34) and (2.36). 

For frozen, equilibrium, and nonequilibrium chemistry, equations 
(2.25) through (2.33), along with appropriate boundary conditions and 
relations for the thermodynamic and transport properties, are the 
governing relations that will be used to describe the viscous shock 


layer. 
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Boundary Conditions 

Conditions at the body surface , -The no-slin boundary conditions are 
used in this study. The surface conditions for n = 0 (see fig. 2) 
are 


u = 0, 


(2.37) 


V = m/p , 


(2.38) 


where the mass injection rate. 


m = (Pv) , 


(2.39a) 


is specified. The mass injection distribution along the body surface 
is assumed to be 


m 


m e 
o 


-3.47 


2 

r 




(2.39b) 


where m^ is the stagnation mass injection rate. This injection 
distribution is the same as that used in reference 9. 

The wall temperature is assumed to be 

T = 
w 


constant . 


(2.40) 
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Also, the surface total enthalpy is given as 


H = 


NS 

Z h.C. . 
, 1 1 


(2. -41) 


Finally, consider the boundary conditions for the elemental and 
species equations (see fig. 2) at the gas-solid interface. The species 
i are transported away from the surface at the rate by diffusion 

and at the rate pvC^ by convection. At the same time, the specie? 
i are being convected to the surface at the rate (pvC^). This flux 
may be considered that from an ablating surface or the injectant from a 
porous solid. The surface species concentrations are given (ref. 20) 
by 


pvC^ + = (^’^^i)-' 


(2.42) 


Since 


pv = (pv)_= m , 


(2.43) 


equation (2.42) may be written as 


(C. - C. .'im = (2.44a) 

V 1 1-/ 1 

and in terms of elemental boundary conditions as 


(2.44b) 
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For the numerical computations, the surface elemental and species 
boundary conditions are written in the following form: 


9w 

3n 


+ Bj^w + $2 


= 0 


(2.45) 


where W equals for the elemental equations and for the 

species equations. For the elemental surface boundary conditions, the 
coefficients and B 2 are defined as follows: 

Multicomponent Diffusion 


p V mn Pr Pr 
o _ s s s s 

Bf 2 _ 

£ y yL 
s 

NS « ,M* NS 9C 

62 ■ -8A- + " a" T 

1=1 1 k=l 

^i 


(2.46a) 


(2.46b) 


Binary Diffusion 


p V nm 

81 = - -a- 

s P3MC 


(2 . 46c) 


B 2 


^1^^£- 


(2.46d) 


Accordingly, for the species boundary conditions, the coefficients Bj^ 
and B 2 are defined as follows: 
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Multicomponent Diffusion 


^1 = 


^2 = 


p V mn Pr Pr 
s s s s 

NS 9C 

- 6,C. + ~ T. Ab., i~ 

1 i- L. , I ik 3n 
1 k=l 

H 


(2 Ala) 


(2.47b) 


Binary Diffusion 


p V ™ Pr Pr 
o _ s s s s 

2 - 
e 


(2.47c) 


62 = - 3 iC,_ 


(2.47d) 


When the diffusion is binary and no injection occurs, the elemental 
concentrations are constant and equal to the free-stream value. For 
this case, the elemental equations are superfluous and the surface 
conditions become 



(2.48) 


However, when the diffusion is multicomponent the elemental distribution 
along the normal to the surface is not constant and the element conser- 
vation equations must be used to calculate the elemental distribution. 
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Conditions at the Shock . -The conditions imposed at the shock were 
calculated by using the Rankine-Hugoniot relations where the flow was 
assumed to be in chemical equilibrium or frozen at the free-stream 
composition. The nondimens ional shock relations are as follows: 


p v" = - sin a , 
•^s s 


(2.49a) 



cos a , 


P = 
s 




P 

s 


h 

s 


p T R* 
s s 

M* C* 


NS 

L h.C. 
i=l " " 


(2.49b) 

(2.49c) 

* 

(2.49(1) 

(2.49e) 

(2.49f) 


and 



(2.49g) 


Also, for the calculations where the chemistry was frozen across the 
shock, 


C. 

r 


s 



(2.49h) 
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When the flow field chemistry is equilibrium, additional relations are 

required to relate the species compositions to the state variables. 

This is accomplished by using a free-energy minimization equilibrium 

program that will be discussed later. 

The transformations used to express the shock— oriented velocities 

u" and v" in terms of the body-oriented coordinate system (fig. 1) 
s s 

are 

u = u" sin(a + 3) + v" cos (a + 3) (2.50) 

s s s 

and 

V = -u", cos (a + 3) + v" sin (a + 3). (2.51) 

s s s 

The normalized shock conditions become 

u=T=H=v=p=P=1 (2.52) 

at n = 1. Note that the elemental and species mass fractions are 
not normalized with their respective shock values. 

Surface Transport 

The rate at which heat is transferred to the surface boundary, q, 
is given by the sum of the convective, conductive, and diffusion contri- 
butions. This is true if radiation and surface reactions such as 
sublimation, vaporization, oxidation, and melting are neglected, as was 
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the case for this study. The net nondimens ional heat transfer to the 
solid interior is given by 



The Stanton number is given by 


St 



00 w 


( 2 . 54 ) 


and the skin friction coefficient is given by 


C 


f 


2e" 



( 2 . 55 ) 
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III. THERMODYNAMIC AND TRANSPORT PROPERTIES 


The thermodynamic properties and h^ and the transport pro- 
perties K^, and £f required for each species considered. 

Also, for the equilibrium calculations, the free energies, F^, are 
required. Since the multicomponent gas mixtures are considered to be 
mixtures of thermally perfect gases, the thermodynamic and transport 
properties for each species are calculated by using the local tempera- 
ture and pressure. Then the properties for the gas mixtures are 
determined in terms of the individual species properties. The methods 
used to calculate the thermodynamic and transport properties are given 
in this chapter. Since all the expressions presented are in terms of 
dimensional quantities, the starred superscript v:ill not be used to 
denote a dimensional quantity. 


Thermodynamic Properties 

Data for thermodynamic properties were obtained from references 23 
and 24. Reference 23 presents thermodynamic properties for 210 
substances for temperatures up to 6000°K. Reference 24 presents 
thermodynamic properties for a smaller number of species than reference 
23 but for temperatures up to 15000°K. The specific heats are ' 
expressed in terms of temperature according to the following polynominal 
equation: 



a, + a-T + a_T^ + a.T^ + a^T^ 
1 2 3 4 5 


R 


(3.1) 
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The coefficients through a^ are given in table I for each 


species considered in this study. In addition to specific heat data, 
the following quantities are 

J- 

given in table I. The enthalpy for species i is obtained from the 
relation 


:udy. In addition to specific heat d^ 
required: 

Values for these quantities are also 


h. = C dT + h^. 

i 298 p. 1 


(3.2) 


298 


where h°. is the enthalpy of formation. Development of the 

1298 

expressions for the free energies, F^, will be given in Chapter IV. 
The mixture frozen specific heat and the mixture enthalpy are given by 
equations (2.13) and (2.15), respectively. 


Transport Properties 

The gas mixtures considered in this work were assumed to be mixtures 
of thermally perfect gases. Accordingly, the expressions for calcu- 
lating the transport properties of the pure species will be presented 
first, followed by the expressions for calculating the transport pro- 
perties of the mixture. 

To a first approximation (first-order kinetic theory), the viscosity 
for a perfect gas (ref. 22) is given as 


2.6693X10 


^ = 




1 


(3.3) 
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and the thermal conductivity for a monotonic gas is given as 


1.9891X10 


-4 


K. 




mono 


1 


or 


K. 

1 


mono 




(3.4a) 

(3.4b) 


whereas the thermal conductivity for a polyatomic gas has an additional 
contribution due to the transfer of energy between translational and 
internal degrees of freedom. Consequently, the expression for the 
thermal conductivity of polyatomic molecules is given (ref. 25) as 



(3.4c) 


In the foregoing equations 

T = temperature in °K, 

= molecular weight, 

a. = collision diameter in X 
1 

= viscosity in g/cm sec, 

= thermal conductivity in cal/sec-°K, 

and 

(2 2 ) 

’ = collision integral for viscosity. 

Expressions for the multicomponent viscosity and thermal Conductivity 
according to rigorous kinetic theory are given in reference 22, However, 
these expressions are cumbersome and computatlonly time consuming. 
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Consequently most studies use approximate techniques for evaluating the 
mixture viscosity and thermal conductivity. 

The mixture viscosity is obtained by using the semiempirical formula 
of Wilke (from ref. 18) 


where 


U 


NS 

E 

i=l 




NS 

j=l 



(3.5a) 


(3.5b) 


The mixture thermal conductivity is obtained by a method analogous 
to that used for viscosity. The relation used (ref. 18) for calculating 
the mixture conductivity is 


K = 


NS 

E 

i=l 


Vi 


NS 

E x.$. . 
j=l J 


(3.6) 


where the $ . coefficients are identical with those that appear in 
the viscosity equation (eq. (3.5b)). 

The binary diffusion coefficients are obtained by using the relation 
given in reference 22 
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2.628X10 T 




V 


3^3/2^/ i 3 


IJ 


2M.M. 
LJl 


(3.7) 


,(1,1) 


is the collision integral for diffusion. The collision 

(1,1) o(2.2) 


where 9, 

integrals and 

reduced temperature. 


are functions of the non dimensional 


= 


e . . /k ’ 
13 


(3.8) 


where is the maximum energy of attraction between colliding 

molecules and k is Boltzmann's constant. The parameters and 

are estimated satisfactorily by combining as 


a. . = ^ [a . + a .1 
IJ 2 L 1 jJ 


(3.9) 


and 


The force constants and e^/k for each species considere3'~ln-'~- 

this study are given in table II. The force constants are those given 
in references 25 and 26. The collision integrals were obtained from 
reference 22, pages 1126-1127, where they are tabulated as a function 
of reduced temperature t"^. These collision integrals are based on the 
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Lennard Jones (12-6) potential for the interaction of colliding mole- 
cules . 

The multicomponent diffusion coefficients are dependent on the con- 
centration of the species whereas the binary coefficients are concentra- 
tion independent. The multicomponent coefficients are obtained by using 
the relations given in reference 1, which may be written as: 


I M. - 

D. . = K. . - jr- K. . 
ij iJ 11 


(3.11a) 


where the quantities K. . are coefficients in a matrix which is the 

iJ 

inverse of the matrix with the following coefficients: 




M. 

3 


NS 

Z 

Jl=l 


for i 7 ^ j 


for i = j 


(3.11b) 


For the binary calculations, it is assumed that all species have 
diffusion coefficients that are equal to an effective binary coefficient. 
The effective binary coefficients are determined for a pre-assigned 
diffusing pair of species. For example, the effective binary diffusion 
coefficients for air will be that of atomic oxygen diffusing into 
molecular nitrogen evaluated at the local temperature and pressure 
according to equation (3.7). However, for the calculations with multi- 
component diffusion, the diffusion coefficients, , are evaluated 

according to equation (3.11). 
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IV. CHEMICAL REACTION MODELS 

Analyses of chemically reacting flows are frequently simplified by 
considering the limiting cases of frozen and equilibrium chemical 
behavior. However, the question naturally arises: Which of these two 

limiting cases is more descriptive of reality? Of course, to answer this 
question the detailed finite-rate behavior of the chemical reactions must 
be accounted for. Yet, an analysis of a complex chemically reacting gas 
mixture, where the chemistry is nonequilibrium, presents problems that 
are not associated with either frozen or equilibrium treatments. These 
problems are primarily those of defining a realistic reactions model, 
of obtaining the appropriate rate constants, and of obtaining a numerical' 
solution. In general, it is the latter problem that precludes a non- 
equilibrium treatment of reacting flows. 

In this study, the finite- rate treatment of the chemical reactions 
will be applied only to a system of neutral air species (0, 0 ^, N, N 2 , 
and NO). However, the equilibrium treatment includes neutral air 
species as well as the complex chemical systems that result from 
injection of water or ablation products into a reacting air stream. 

The methods used to make the equilibrium and finite rate calculations 
will presently be discussed. 

Equilibrium Chemistry 

The mixture equilibrium composition is determined for a given 
temperature and pressure by the free-energy minimization technique. 
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Reference 27 describes the free-energy minimization computer program 
used in this study. 

In the free-energy minimization method, the equilibrium constants 
are not used. In fact, specification of the various reactions is not 
required. However, the species for which the equilibrium composition 
is to be determined must be specified and the elemental composition 
must be known. Also, the free energy for each species is calculated 
from which the total free energy of the mixture is determined. This 
total free energy is then made a minimum with respect to any possible 
change in composition for the given temperature and pressure. Since the 
criterion for equilibrium at constant temperature and pressure is that 
the change In free energy be zero, dF = 0, the equilibrium composi- 
tion is determined when the total free energy of the mixture is a 
minimum. 

The free energy for a mixture of thermally perfect gases is obtained 
by using the following definition for free energy of species i: . 

F. = h. - TS. , (A.l) 

111 

where S is the entropv of species i. The standard state (pure 
1 

substance at one atmosphere of pressure) free energies. 



(4.2) 
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are a function of temperature only. The enthalpies are obtained 
relative to a base state, taken as the elements in their most natural 
form at 298°K and one atmosphere (ref. 23). For a constant pressure, 
an expression for the entropies in terms of specific heats is obtained 
from the combined first and second laws of thermodynamics as 


C dT 



(4.3) 


or 


^ i = ^ i + ^ 
298 


C dT 

Pj 
298 T 


(4.4) 


With 


h°. = C dT + h". 

X 298 p. 1^98 


(4.5) 


equation (4.2) may be written as 


= F°. - - F°. 1 

^ ^298 1_ Jl ^298 ^298j 

C dT 

-^298 ^ ^298 ~T ‘ 

X 


(4.6) 


For constant temperature, the change in free energy of a thermally 
perfect gas is given by the expression 
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dF^ = RTdUnp^ (-^.7) 

By integrating equation (4.7) at constant temperature from a standard 
state where the free energy and pressure are and one atmosphere, 

respectively, to some other state F^ and p^ , the following is 
obtained: 


or 


since 


F. = F°. + RT£np. (4.8a) 

11 ^1 



(4.8b) 


p. = px. (4.9) 

^1 1 

and X denotes the mole fractions. The expression for the free 
i 

energy of the mixture becomes 

NS ■ , 

F = E x^ |^F°^ + RT (Jlnp + S^nxJJ (4.10) 


where the free energies in the standard state are given by equation (4.6). 



41 


By minimizing the above equation subject to the constraint equations 


NS 

E 

i=l 





M 


X. 

1 


(4.11) 


the equilibrium composition is obtained. Details of the solution 
procedure are given in reference 27. 


Finite Rate Chemistry 

When chemical reactions proceed at a finite rate, the rate of pro- 
duction terms, are required. The production terms appear in the 

energy equation (eq.(2.4)) when formulated in terms of temperature and 
in the species continuity equations (eq. (2.5)), For a multicomponent 
gas with NS chemical species and NR chemical reactions, the chemical 
equation describing the overall change from reactants to products may 
be written in the general form 


NJ 

Z 

i=l 



a. X 

i 



NJ 

Z 

i=l 


3. X. 
i.r 1 



(4.12) 


where r = 1, 2, . . .NR and NJ is equal to the sum of the species and 

the catalytic bodies. The quantities a. and 3. are the stoich- 

1 , r a j t 

iometric coefficients for reactants and products respectively whereas 

and k, are the forward and backward rate constants. Also, 
f,r l),r 

the quantities X^ denote the concentrations in moles per volume of 
the distinct chemical species NS and the catalytic bodies. The 
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catalytic bodies (NJ-NS) may be chemical species or linear combinations 
of species that do not undergo a chemical change during the reaction. 
The rate of change of any species as a result of a particular reaction 
(ref. 28) is 



i.r 


- a. 



NJ 

n 

j=i 




NJ 

n 

j=i 


X. 

j 



(4.13) 


The above equation may be rewritten as 


/ \ r r NJ a. NJ . 3 . n 

[«i,r - (^/) - A.r (’'/I 


where the mole-mass ratio, Yj . is defined as 


X. C. 

Y = -J. = -J- 
P. M. 
1 J 


(4.15) 


and 


w . 

1 


^dXi 

M. dt 


dPi 

dt 


(4.16) 


In order to find the net mass rate of production of the i species 
per unit volume, equation (4.14) must be summed over all reactions r. 
Thus, one obtains 


w. = M. 

1 1 



NJ 

n 

j=i 



ot. 

J.r 
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which is the general rate equation for a gas mixture. 
The mole-mass ratios for the NS species are 


y. = c./M. 

J 1 2 


j = 1, 2, . . . NS, 


(4.18a) 


whereas for the catalytic bodies, the following expressions (ref. 2) 
are used: 




NS 

= E Z 




j = NS + 1, 


NJ. 


(4.18b) 


The constants ^ are determined from the linear dependence of 

the catalytic bodies upon the NS species. Values of these constants 
are given in reference 11 and in table III. 

The chemical reactions used in the study are as follows: 


r = 1 0„ + M, r ^ 20 + M, (4.19a) 

2 1 <«i 1 

2 N 2 + ^ 2N + (4.19b) 

3 N^ + N i 2N + N (4.19c) 

4 NO + > N + 0 + (4.19d) 

5 NO + 0 .2 O 2 + N (4.19e) 


N 2 + 0 ^ NO + N . 


6 


(4.19f) 
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The reaction constants for these equations are expressed in the modi- 
fied Arrhenius form, where the forward rate is given as 


C2 


k. = T 
f 


^ exp j^; 


ilnCO^ - Clj. X 10 /tJ 


and the backward rate is given as 


r 3 1 

k^ ^ = T exp UnDO^ - Dl^ X 10 /TJ 


(4.20a) 


(4.20b) 


The experimentally determined coefficients for these equations are those 
given in reference 29 and are tabulated in table IV. Therefore, for a 
specified temperature, density, and species composition, equation (4.19) 
along with equations (4.18), (4.19), and (4.20) are used to determine 
the rates of production for a multicomponent gas. The manner in which 
the production terms are written for numerical solution was discussed 
in Chapter II. (See equations (2.34) and (2.36)). 
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V. METHOD OF SOLUTION 

In this chapter the procedure for solving the viscous shock-layer 
equations is presented. Development of the stagnation streamline 
equations is presented first where the nonlinear set of partial differ- 
ential equations is reduced to ordinary differential equations by using 
a truncated series method. Then the finite-difference expressions 
used to transform the differential expressions to algebraic expressions 
are presented. The solution procedure is then discussed. This procedure 
is essentially that developed by Blottner and Flugge-Lotz (ref. 30) 
for compressible boundary layer solutions and applied by Davis (refs. 

7 and 8) to the viscous shock- layer problem. 

Basically, the solution procedure is to assume an initial stagnation 
profile for all dependent variables, do iterative calculations locally 
until the solution converges, advance to the next body station using 
the previous station profiles as an initial guess, and repeat the 
procedure until around the body solution is obtained. Once this is 
accomplished, the solution may be iterated in ^ to remove the 
assumptions (thin shock- layer form of the n-momentum equation and 
assumptions regarding the shock geometry) that were used in the first 
approximation . 

Stagnation Streamline Equations 

When downstream numerical solutions are required, it is necessary 
to have an accurate solution for the flo^'/ along the stagnation 
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streamline. A truncated series, which has the same form as that used 
by Kao in reference 31, is used to develop the stagnation streamline 
equations. The flow variables are expanded about the axis of symmetry 
with respect to the nondimensional distance ^ along the body as 
follows : 


P(C,n) 

= P^(n) + + — , 

(5.1a) 

u(?,n) 

= u^(n)C + — , 

(5.1b) 

v(?,n) 

= v,(n) + — , 

(5.1c) 

p(C,n) 

= pj^(n) + — , 

(5. Id) 

T(?,ri) 

= T^(n) + — , 

(5.1e) 

h(^,n) 

= h^(n) + — , 

(5. If) 

p(5,n) 

= + — , 

(5.1g) 

K(?,n) 

= K^(n) + — , 

(5.1h) 


(5.1i) 

(5.1j) 


c (^,n) = c (n), + — 

p ^1 


(n) + — 

j 
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and 

c.(C,n) = c (n) + — . (5.1k) 

1 

The shock stand-off distance is written as 

"s ■ "is "2s^^ + — • (5.2) 

Furthermore, ^ is small and the curvature , ic , is approximately 
one in the stagnation region. Consequently, the geometric relations 
(see figure 1) including terms of order ^ can be written as 

6-5 (5.3) 


and 



Therefore , 


sin(a + B) ~ 1 


(5.4) 


(5.5) 


and 


cos (a + B) ~ 




1+n 


Is 


(5.6) 
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The shock relations (equations (2.49) in terms of expanded variables ) 
become 




— 


, _ 1 

2 \ 

ym \ 

‘ oo 

1 

"D j 

ro 1 

E-t) ( 


1 

. l^ls ) 

J ■ 


and 



(5.7d) 


along with equations (2.49f) through (2.49h). 

By examining the above equations it is seen that the equations for 

u and P contains n„ . This term cannot be determined from the 
s s 2s 

stagnation solutions since it is a function of the downstream flow. 

Consequently, a value must be assumed for i'2s’ this study, it 

was assumed to be zero to start the solution but this assumption is 

then removed by iterating on the solution by using the previous shock 

stand-off distances to define n„ . The effect of the downstream shock 

2s 

shape on the stagnation point solution is elliptic rather than parabolic. 

By substituting the above relations into the governing viscous 
shock-layer equations, the following relations are obtained: 
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d^w . dw , 

TT “i cFi V “3 = 0 
dn 


(5.8) 


where w assumes the values of u^^, and for the 

s-momentum, energy (enthalpy and temperature), species continuity, and 
elemental continuity equations, respectively. The coefficients are 
defined as 

s-momentum 


1 '*''1 , “is "is^ls^'ls '^l’'l 

“i ■ - *r ^ ^ 1 : 5 ^ ■ , 2 „ - If 

“l “ “is “1 


(5.9a) 


ct„ = - 


*^is r ^1 


'^Is ^Is’^ls'^ls '^1^1 

+ . . , (J + 1 ) + 


1+r.i^n 


E^U 


Is 


^Is^ls^ls 

“^“is l^'""ls“l “1 


p] 


a„ = 


" 2P, n , 
Is Is 


l^'^ls'’! “ls“ 


i 


P P n^ n dP 

P + 2s 1 _ _ls 

2 P 


Is 


n^ dn 
Is 


dn J 


(5.9b) 


(5.9c) 


Energy (temperature) 


1 ‘**^1 , I n “is 
“1 ■ iT *T 

1 Is 


n, NS 

Z J C 

K, K, i=l ^1 ^i. 
Is 1 1 


n- p, C V pc V 
Is Is p^g s 1 p^ 1 


e^K 


Is 


K. 


(5.10a) 
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= - 


n 1 w „w„ 
Is s2 2 

Is 1 


a , = 


'"l 

e ^ T , K K 
Is Is 1 


2 . 

^ Is'^sl T 
^ ^ w 

e x, K, K, 

Is Is 1 


Energy (enthalpy) 


_ ^ dPi ^ dPr^ 

■ n/’' ■ Pr ^ “1 " 


n, p, Pr, V- p,Pr,v, 
Is^^ls Is Is 11 


2 

e y 


Is 


a-B = 0 


= 


'■' islia !!! Q 

i. H, - n. 

Is Is y^ [_ 1 


dt ^ ( H - l )' l ' 
Is ^^''-'is 




Species Continuity 


a, = 


1 dPL 


1 PL dn 


+ 


b 


Pls'ls ^lsPj^l 

e^PL 


Oi ^ = - 


2 .1 
n T w . 
Is X 

2 

e PL 


= 


k [: 


dPH ,,,,,, '■^Is , 


i ] 


and 


Elemental Continuity 


(5.10b) 

(5.10c) 

(5.11a) 

(5.11b) 

(5.11c) 

(5.12a) 

(5.12b) 

(5.12c) 
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1 dPL . "is "is^ls^'ls^l'^l 


e^PL 


(5.13a) 


a2 = 0 


(5.13b) 


a- = 

3 ^ 

PL 


a 


PM 

+ (j + 1) 


dn 


1+nn 


Is 


] 


(5.13c) 


The remaining equations may be written as follows: 
Global Continuity 


dn "is^) ^ ^Is^ls^l^l]" 


- (j + l)n^^ (l + n^^n) ^ PisU^sPiUi 


n-momentum 


- P-.V. 


dn 


P, dn 

Is 


(5.14a) 


(5.15a) 


When the thin shock- layer approximation is made, the n-momentum equation 
becomes 



(5.15b) 


Also, the ?2 term that appears in equation (5.9c) can be expressed in 
terms of "one" quantities by using the follovring equation: 
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AW ^ + B W + C W = D (5.16) 

n m, n-1 n m,n n in, n+1 n 

The subscript n denotes the grid points along a line normal to the 
body surface, whereas the subscript m denotes the grid stations along 
the body surface. Equations (5.16) along with the boundary conditions 
constitute a system of the tridiagonal form, for which efficient 
computational procedures are available (ref. 32). 
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Solution Procedure 

Once the free stream conditions and body geometry are specified, 
initial profile values for the dependent variables are assumed along 
with a guess value for the stagnation shock stand-off distance. Then 
the shock quantities are calculated by using the shock equations 
previously described. For the equilibrium shock calculations, initial 
guess values for the density and mixture molecular weight are made. 

Then a double iteration scheme is used (see figure 3(a)) to obtain a 
converged solution. For the shock calculations where the chemical 
composition is frozen at the free stream value, only a guess value for 
the nondimensional shock density is made (see figure 3(b)). Then 
iterative calculations are made until convergence is achieved. 

With the shock conditions established, the viscous shock-layer 
equations are solved in the sequence shown in figure 4(a) for the 
equilibrium calculations and the sequence shown in figure 4(b) for the 
frozen and nonequilibrium calculations. Finite difference methods are 
used to integrate the equations. Along the stagnation streamline, 
equations (5.8) through (5.15) are used, whereas for the downstream 
equations (2.25) through (2.33) are used. 

First consider the equilibrium calculations. As shown in figure 
4(a), the elemental continuity equations are integrated to determine 
elemental profiles. Next, the energy equation in terms of total 
enthalpy is integrated to determine a total enthalpy profile. Then 
a temperature profile is calculated iteratively such that is consistent 
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with the enthalpy values resulting from the energy equation. This is 
accomplished by, converting the total enthalpy to static enthalpy values, 
then perfoimiing the following procedure at each grid point across the 
flow field: With local values of pressure, elemental composition, 

and an assumed value of temperature, an equilibrium calculation is made, 
whereby new values for species composition and static enthalpy are 
obtained for the assumed temperature value. Additional iterations are 
made until the temperature corresponding to the local enthalpy is 
determined. Once the temperature profile is calculated, the species 
profiles are also known. 

For the frozen and chemical nonequilibrium calculations, the above 
procedure is somewhat simpler (see figure 4(b)). That is, the species 
continuity equations are integrated to determine species profiles for 
nonequilibrium flow. For frozen flow, the species continuity 
equations are not required. For both nonequilibrium and frozen flow, 
a temperature profile is determined directly by integrating the energy 
equation that is written in terms of temperature. 

The procedure for solving the remaining equations is identical for 
all three chemistry models. With new values for species and tempera- 
ture profiles, a density profile is calculated by using the equation 
of state. Next, thermodynamic and transport properties are calculated. 
Then the s-momentum equation is integrated to determine the tangential 
velocity profile. Thereafter, the continuity equation is integrated 
to determine first the shock stand-off distance and second the v 
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component of velocity. Next, the n-momentum equation is integrated to 
determine the pressure profile, and the equation of state is used once 
again to update the density profile. This procedure is repeated until 
the differences between values of quantities for the current iteration 
and the previous iteration satisfy the convergence criteria. Once the 
solution has converged, the computation advances to the next body 

station using as initial profiles those computed for the previous 

« 

station. Again, iterative computations are made until convergence 
is achieved. Hence, the computation advances in ^ until the desired 
body station is reached. 

The foregoing procedure provides a first approximation to the flow 
field solution since the following assumptions are used: 

a. The thin shock-layer form of the n-momentum equation (2.32b) 
was used. 

b. The stagnation streamline solution was independent of downstream 

influence (approximation of local similarity where n„ = 0) . 
dn ' ® 

g 

c. The term was equated to zero at each body station. 

d. The shock angle a was assumed to be the same as the body 
angle 0. 

These assumptions are then removed by making a second approximation to 
the flow field solution whereby the v components of velocity from 
the first approximation are retained and averaged with the local and 
current v velocities and then used in the n-momentum equation (eq. 

2,32a) for determining the pressure profile. Furthermore, the shock 
stand-off values from the first approximation are used to evaluate dn^/d? 
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which is then used in the subsequent approximation. If necessary, 
additional iterations are performed until the total solution converges. 
Hence, the viscous shock- layer equations are solved as parabolic equa- 
tions, yet retaining effects, the nature of which are elliptic and 
hyperbolic. 
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V. RESULTS AND DISCUSSION 

• ^ 

Introduction 

In this chapter numerical solutions to the previously discussed 
viscous shock-layer equations are presented and discussed. The solu- 
tions are those for flow about a 45° hyperboloid at zero-degree angle 
of attack. First, comparisons of present results are made with data in 
the literature. Second, results are presented which demonstrate the 
effects of chemistry model, diffusion model, surface cataly ticity , and 
mass injection on flow parameters and surface transport. 

All computations, unless otherwise indicated, were made by using 
a constant step size of 0.2 in the ^ -direction and a variable step 
size in the n -direction where a total of fifty grid points were 
used. Also, two iterations were made in the ^ -direction. Table V 
summarizes the altitude, velocity, and free stream properties used. 

The free stream properties are those given by the U. S. Standard 
Atmosphere (ref. 33). 

All computations were made by using a Control Data Corporation (CDC) 
6600 computer. The program language was FORTRAN IV. The maximum 
storage requirement of the different programs was 147000 core locations, 
and the program was dimensioned to accommodate fifty grid points along 
the normal to the body, fifty grid points along the body surface, four 
chemical elements, and twenty chemical species. A plotting routine 
was also included in the program. Table VI summarizes the total time 
required to calculate various cases and the number of iterations 
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required at the stagnation station and a representative downstream 
station. For purposes of comparison, time and iterations required for 
an ideal gas calculation are also included. 

For each body station, the convergence criteria was that the relative 
difference be less than 0.001 for both the temperature and tangential 
velocity derivatives at the wall. This criteria was used for all 
calculations except for those where the temperature gradient was 
approaching zero due to mass injection. In addition, a minimum of 
four iterations was required at each body station. 

The equilibrium flow solutions require large computing times because 
of the time required to make a chemical equilibrium calculation. For 
equilibrium flow calculations where the diffusion is binary, a major 
reduction in computing time can be realized by using tabulated or 
curve fitted data for the chemical composition. This is also appli- 
cable for mass Injection provided the injectant has the same elemental 
composition as the free stream. However, this approach was not used 
in the present study since most of the equilibrium calculations were 
for multicomponent diffusion or binary diffusion with an injectant 
whose elemental composition was not the same as the free stream 
composition. 

Comparison With Previous Analyses 
In order to assess the results obtained for the present analyses, 
comparisons were made of the results for each of the three chemistry 
models with data reported in the literature. No attempt was made to 
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obtain exact comparisons since the solution is dependent upon step size, 
values of transport and thermodynamic quantities, details of the given 
chemistry model, etc. 

Frozen and equilibrium calculations were compared with the data of 
Whitehead as reported in ref. 9. The comparisons were for the following 
conditions: a free stream velocity of 20,000 ft/sec, an altitude of 

200,000 ft in the earth’s atmosphere, a 45° hyperboloid with a one- 
inch nose radius, and a surface temperature of 1800°R. Even though 
the solution procedure used herein is essentially the same as used by 
Wliitehead, several differences exist. Whitehead results accounted for 
shock slip at the outer boundary condition and slip and temperature 
jump at the surface, whereas the present analyses use the Rankine- 
Hugoniot relations to describe the shock conditions and assumes no 
slip or temperature jump at the wall. Furthermore, the expressions used 
to calculate the thermodynamic and transport properties are not the 
same as used by Whitehead. 

Figure 5(a) through (c) show comparisons of the present analysis 
with that of Whitehead's for shock stand-off, skin friction coefficient, 
and nondimensional heat transfer rate where the chemistry is frozen 
at the free stream composition. Table VII presents a detailed listing 
of the shock and wall values that result from the present calculation. 
The heat transfer results are in excellent agreement; however, the 
agreement is not as good for shock stand-off distances or skin friction 


coefficients . 
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Figures 6(a) through (c) show the same comparisons for equilibrium 
air where the Lewis number is assumed to be one. Table VIII presents 
a detailed listing of the stagnation shock and wall values for the pre- 
sent calculation. Whitehead's equilibrium results were obtained 
assuming air to be represented as one "effective" species with a 
constant molecular weight equal to the free stream value whereas the 
present analysis treated air as a multicomponent mixture of five species. 
The present calculation predicts greater values for shock stand-off 
distance (as would be expected due to Whitehead's assumption concerning 
the molecular weight), skin friction coefficient, and heat transfer. 

Results of the equilibrium air calculations were also compared with 
the data of Edelman and Hoffman (ref. lA) V7hich was for the stagnation 
station only. Comparisons of stagnation temperature and species 
concentration profiles are presented in figures 7(a) through (d). 

Figures 7(a) and (b) show a comparison of the temperature and species 
profiles, respectively, for a 280,000-ft altitude, a 26,000-f t/sec 
free stream velocity, a 12— inch nose radius, a 2A30°R wall temperature, 
and a Lewis number of 1.4. Table IX presents a detailed listing of 
the stagnation shock and wall values for the present calculation. 

Figure 7(c) and (d) show the same comparisons for a 150, 000-ft altitude, 
a 22,520-ft/sec free stream velocity, a 2657°R wall temperature, a 
0.5- inch nose radius, and a Lewis number of 1.4, Table X presents the 
stagnation shock and wall values corresponding to the above conditions 


for the present calculation. 
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When the results of the two calculations for both sets of conditions 
are compared, two differences are apparent: first, the shape of the 

temperature profiles, and second, the shock-layer thickness. The 
curvature of the temperature profiles as calculated by reference 14 
decreases monotonically with increasing distance from the surface. 
Temperature profiles res.ulting from the present equilibrium calculations 
are not smooth, especially at the higher altitude (lower pressure) 
conditions. However, it is believed that the temperature profiles 
resulting from the present calculations can be substantiated by 
considering the chemical reactions that occur within the shock layer. 

A more detailed discussion of this assessment is given in the next 
section of this chapter. 

The differences in shock layer thickness are expected since the 
results presented in reference 14 were for the stagnation streamline 
only and, therefore, did not account for the downstream influence on 
shock stand-off distance. However, when the dc^^^nstream influence on 
shock shape is accounted for as in the present analysis, a larger 
shock layer thickness results. This effect was shown in reference 7 for 
an ideal gas . 

Results of the nonequilibrium calculation were compared with the 
stagnation viscous shock-layer results of Blottner (ref. 11). Compari- 
sons were made of temperature, tangential velocity, and species concen- 
tration profiles. 

For the chemical kinetics model, Blottner considered two additional 
species (NO"*" and e~) and one additional reaction equation (N+OtNO*^ + e") 
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however, for the remaining reactions, identical rate expressions 
were used. The calculations are for a 200,000 ft altitude, a 20,000 
ft/sec free stream velocity, a 1800°R wall temperature, and a catalytic 
wall. Both analyses use the Rankine-Hugoniot relations to describe the 
shock conditions and use the no-slip or temnerature jump conditions at 
the wall. Also, for purposes of comparison, the present calculation 
was made with a constant Lewis number of 1.4. 

For two reasons, differences in shock— layer thickness are expected. 
First, the results of reference 11 do not account for the downstream 
influence on the stagnation shock stand-off distance as does the present 
analysis. When the downstream effect is neglected, a significant under- 
prediction of stagnation shock stand-off distance can occur (refs, 7 
and 34). Second, the analysis of reference 11 did not include the 
effect of the shock-layer thickness on the shock radius of curvature 
as did the present analysis. Reference 34 (p. 91) found that when the 
shock-layer thickness contribution to the radius of curvature is 
neglected, an underprediction of the stagnation shock stand-off distance 
occurs. 

Table XI presents a detailed listing of stagnation shock and wall 
conditions for the present calculation (a constant Lewis number of 1.4). 
Figure 8(a) shows a comparison of temperature profiles. The most 
apparent difference is in the shock stand-off distance. Nevertheless, 
near the wall the temperature profiles are in good agreement. The same 
can also be said for the tangential velocity profiles as evidenced 
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by figure 8(b) where the agreement is good near the wall. 

Figure 8(c) shows a comparison of the species concentration profiles 
for the two analyses. Appreciable differences exist. However, it is 
believed that these differences are due primarily to differences in 
shock stand-off distance. The shock stand-off distance for the present 
analysis is approximately 22 percent greater than that of reference 11. y 

Consequently a larger distance is available for reactions to occur 
in particular the dissociation reactions. This is evident in figure 
8(c) which shows that, for a given value of n, the amount of dissocia- 
ted flow predicted by the present calculation is greater than that 
predicted by reference 11. , 

Chemistry Models 

In this section results of calculations for a free-stream velocity 
of 20,000 ft/sec, an altitude of 200,000 ft in the earth's atmosphere, 
a one- inch nose radius, and a surface temperature of 2700 R using 
frozen, equilibrium, and nonequilibrium chemistry are presented. Multi- 
component diffusion, no mass injection, and a noncatalytic wall are 
assumed. Surface catalyticity effects are discussed in the next 
section. First, velocity, temperature, and species concentration 
profiles are presented for each chemistry model. Next, the equilibrium 
temperature profiles are shown and their dependence on chemical 
reactions and pressure variations are discussed. Then the effect of 
each chemistry model on tangential velocity, temperature, and chemical 
species profiles, heat transfer, Stanton number, skin friction, shock 
stand-off distance, and wall pressure is shown. 
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Figures 9 through 12 show the velocitv, temperature, and species 
profiles at various body stations for each chemistry model. There are 
several features of these profiles that are of interest. The tangential 
velocity profiles are practically linear in the outer region of the 
shock layer, that is, the flow has an outer region of essentially 
uniform vorticity. This is in marked contrast to the classical velo- 
city boundary- layer profile that is derived with a zero velocity 
gradient imposed at the outer boundary. However, the flow considered 
here is in the viscous- layer regime (ref. 35) where the viscous effects 
extend over an appreciable fraction of the shock layer, thereby invalida- 
ting the boundary-layer concept. 

The equilibrium temperature profiles are shown in figure 10(b). 

Near the wall, the equilibrium temperature profiles do not show the 
monotonically decreasing curvature associated with frozen and non- 
equilibrium profiles. The equilibrium temperature profiles have two 
distinct Inflection points. From the body to the shock, the profiles 
are characterized by concave, convex, and concave segments. The 
curvature of the convex segment diminishes as the flow moves downstream. 
To the author's knowledge, this complex shape of the equilibrium air 
temperature profiles has not been observed previously. 

The bulges appearing in the temperature profiles can be explained 
in terms of chemical reactions. This is demonstrated in figure 13, 
where stagnation temperature profiles and species concentration profiles 
are shown. Temperature profiles for an equilibrium chemistry solution 
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and a frozen chemistry solution (where the chemistry is frozen at the 
equilibrium shock composition) are shown. The temperature profile for 
the equilibrium calculation is much fuller than the profile for the 
frozen calculation. This should be, of course, because exothermic 
reactions occur near the wall. Essentially all chemical reactions occur 
within that half of the shock layer which is adjacent to the wall; that 
is, essentially all the atom recombinations occur near the wall. The 
recombinations also occur over two distinct and separate regions. The 
oxygen recombines between the wall and r) = 0.11. This causes the 
inner bulge in the temperature profile. Then there is a small region 
(r| of 0.11 to 0.16) where the recombination process is negligible. 
Finally, practically all the atomic nitrogen recombines between n 
values of 0.16 and 0.5. This causes the outer bulge in the temperature 
profile. 

The extent to which the equilibrium temperature profiles differ from 
smooth profiles will depend upon the pressure, hence the altitude for 
a given free stream velocity. This is demonstrated by Hansen in 
figure 1 of reference 36. The data of Hansen show that the recombination 
reactions occur over narrow temperature intervals at low pressures 
(0.0001 to 0.1 atm); however, the temperature intervals required for 
nitrogen and oxygen recombination increase significantly with increas- 
ing pressure. Furthermore, the temperature interval between the 
completion of nitrogen recombination and the onset of oxygen recom- 
bination decreases with increasing pressure. Consequently, as the 
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pressure Increases, the chemical composition experiences a more gradual 
change with temperature. The pressures were 0.1 atmosphere or less 
for the temperature profiles shown in figure 10(b), therefore, recom- 
binations over narrow temperature intervals would be expected. 

To demonstrate the effect of the recombination reactions on the 
equilibrium temperature profiles at higher pressures, an equilibrium 
calculation was made for an altitude of 100,000 ft and a velocity of 
20,000 ft/sec (a detailed listing of the shock and wall values resulting 
from this calculation are presented in table XII.) The stagnation 
pressure was approximately 5.2 atmospheres. Figure 14(a) shows the 
stagnation temperature profile, while figure 14(b) shows the stagnation 
species profiles. From figure 14(a), it is seen that the curvature of 
the equilibrium temperature profile is essentially monotonic, decreasing 
with increasing p. Hence, the results of the higher pressure calcu- 
lation agree qualitatively with the aforementioned effect of pressure 
on equilibrium temperature. 

Figures 11(a) through 11(d) show the nonequilibrium species profiles 
at different body stations for a noncatalytic wall. The chemistry 
across the shock is frozen at the free stream composition (C^ = 0.24; 

= 0.76). The figures show that most of the dissociation occurs 
in the outermost portion of the shock layer. Furthermore, only a 
small amount of recombination occurs near the wall. This is due to 
the low density flow conditions which reduce the effectiveness of 
recombination reactions; consequently, the amount of dissociated flow 
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at the wall is significant (approximately one-fifth of the total mass) . 
This is important because the energy invested in dissociation is not 
recovered at the wall in nonequilibrium flow. 

Figures 12(a) through (c) show the equilibrium concentration profiles 
at different body stations. It is seen that a large percentage of flow 
is dissociated in the outer portion of the shock layer. The amount of 
dissociation decreases rapidly downstream. Also, essentially all the 
dissociated species recombine in a relatively small region near the 
wall with the wall composition being approximately that of the free 
stream composition. This, of course, means that the energy of 
dissociation is recovered. The significance of this in terms of heat 
transfer rate will be discussed later. 

Comparisons of results using the three chemistry models are shown 
in figures 15 through 21. The effect on tangential velocity and 
temperature profiles, shock stand-off distance, heat transfer rate, 
Stanton number, skin friction and wall pressure are shown. Note that 
all the aforementioned quantities are significantly influenced by the 
chemistry model with the exception of the tangential velocity profiles 
(fig. 15) and the wall pressure distribution (fig. 21). Of particular 
significance are the large differences in shock-layer temperatures for 
the different chemistry models. Figures 16(a) and (b) show the 
temperature profiles at an C of 0.0 and 2.0, respectively, for each 
chemistry model. The conditions at the shock were identical for the 
frozen and nonequilibrium calculation since the chemical composition 
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was assumed to be frozen across the shock at the free stream composition. 
For the stagnation streamline (fig. 16(a)), the frozen shock temperature 
is more than twice that of the equilibrium shock temperature. Moreover, 
the equilibrium temperatures across the shock layer are much less than 
those of the frozen or non equilibrium values except in the region near 
the wall. This is due to the large amount of energy invested in 
dissociation rather than static temperature. Also, the nonequilibrium 
temperatures are less than the frozen values because some dissociation 
occurs within the nonequilibrium shock layer. Also at the downstream 
(fig, 16(b)) station, C = 2.0, the temperature differences are 
noticeably less than those at the stagnation station because the 
equilibrixffli and nonequilibrium compositions are more nearly like the 
frozen composition at the lower temperatures. 

Figure 17 shows how the chemistry model influences shock stand-off 
distance. The shock stand-off distances corresponding to the equili- 
brium calculation are considerably less than those for frozen and 
nonequilibrium calculations. Also, the shock stand-off distances for 
the nonequilibrium calculations are somewhat less than those for the 
frozen calculations. These results are expected since the shock- 
layer densities for the nonequilibrium and especially the equilibrium 
calculations are greater than the frozen flow densities. 

Consider next the effect that the chemistry model has on heat 
transfer rate distribution. As is shown in figure 18, the chemistry 
model has a very pronounced effect on nondimensional heat transfer. 

Both the frozen and equilibrium chemistry models significantly 
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over-predict the heat transfer for a noncatalytic wall. This is 
especially so in the stagnation region where the frozen and equilibrium 
values are approximately 50 percent greater than the nonequi librium 
values. This is particularly significant because the entry corridor 
for manned entry vehicles includes the altitude-velocity conditions used 
in the present calculations. Furthermore, the heating rates are 
significant at these conditions, and therefore, the need for non- 
equilibrium analyses is clearly evident from these results. 

When the heat transfer rate is expressed in terms of Stanton number 
(fig. 19), the results for the three chemistry models are in closer 
agreement than were the nondimens ional heat transfer rates (fig. 18). 

As was pointed out previously (fig. 11), a significant amount of energy 
is invested in dissociation at the wall in nonequilibrium flow. Con- 
sequently, the enthalpy at the wall is larger for the nonequilibrium 
and noncatalytic wall calculation and necessarily the enthalpy poten- 
tial is smaller. The maximum difference between the equilibrium and 
nonequilibrium Stanton numbers is about 13 percent and occurs at the 
stagnation point. 

The effect of the chemistry model on skin friction distribution is 
shown in figure 20. This figure shows that the skin friction distri- 
butions for the equilibrium and nonequilibrium calculations are about 
the same. The values for the frozen calculation are greater than 
the nonequilibrium values, but the differences never exceed 9 percent. 

The effect of chemistry model on wall pressure distribution is 
shown in figure 21. Also shown is the Newtonian pressure distribution. 
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The non equilibrium and frozen pressure distributions are identical. 

Also, the equilibrium and Newtonian pressure values are in close 
agreement, yet less than the nonequilibrium pressure. At C = 3.0, 
the nonequilibrium pressure is 10.8 percent greater than the equilibrium 
and the Newtonian pressures. 

Catalytic Wall Effects 

Results presented in the previous section demonstrated how the flow 
field chemistry can influence flow parameters and surface transport. 

This section demonstrates how surface catalyticity influences the same 
flow parameters and surface transport. Either noncatalytic or equili- 
brium catalytic walls are assumed. For the equilibrium catalytic wall, 
the gas composition at the wall is the equilibrium composition for the 
wall temperature and pressure. The elemental composition at the wall 
is governed by multicomponent diffusion and is, therefore, not 
necessarily equal to the free stream elemental composition. 

Comparison of the noncatalytic and equilibrium catalytic wall 
calculations is made for a 2700°R wall temperature, a 20,000 ft/sec 
free stream velocity, and a 200,000 ft altitude. The diffusion for 
both calculations is multicomponent. Table XI provides a detailed 
listing of the stagnation shock and wall values for the noncatalytic 
and equilibrium catalytic calculations. 

The wall catalyticity has negligible effect on temperature and 
velocity profiles, and therefore, comparisons are not presented. 

However, the species profiles for the flow configuration examined were 
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necessarily influenced by the wall catalyticity. This is shown in 
figures 22(a) through (c) where comparisons are made at an C of 0, 

1.0, and 3.0, respectively, of nonequilibrium species profiles with 
noncatalytic and equilibrium catalytic wall conditions. It is seen 
that the species profiles for the Innermost 70 percent of the shock 
layer is influenced by the wall catalyticity. 

Figures 23 through 26 show how wall catalyticity influences shock 
stand-off distance, heat transfer, Stanton number, and skin friction. 

With the exception of heat transfer, the wall catalyticity has little 
effect. The maximum relative difference between the catalytic and non- 
catalytic calculations were: 3,5 percent for shock stand-off, A. 9 

percent for skin friction, and 9.2 percent for Stanton number. The 
maximum relative difference in heat transfer was 48 percent. 

For purposes of comparison, figure 24 also Includes the heat transfer 
distribution where equilibrium flow field chemistry was used. It is 
seen that the nonequilibrium equilibrium catalytic wall results and 
the equilibrium noncatalytic wall results are in close agreement. 

For both of these calculations, the energy invested in dissociation is 
recovered. However, in the nonequilibrium noncatalytic wall calculation, 
an appreciable amount of dissociation is present at the wall (see 
fig. 11), and, hence, the wall heat transfer is substantially reduced. 

Diffusion Models 

Prior to this study, stagnation and downstream viscous shock- layer 
results that treated the diffusion as multicomponent were not available. 
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In this section, comparisons are shown between multicomponent and binary 
diffusion models. The effects on flow parameters and surface transport 
are presented for both equilibrium and nonequilibrium air. 

The binary approximation, as indicated previously, is that all 
species have the same diffusion coefficient, which is assumed to be the 
coefficient for molecular nitrogen diffusing into atomic oxygen. This 
approach is used to evaluate the binary diffusion coefficient because 
atomic oxygen and molecular nitrogen are the dominant species in the 
shock layer. (See figures 11 and 12.) The multicomponent diffusion 
coefficients are determined according to equation (3.11), 

The flow conditions were as Indicated previously that is, a 200,000- 
ft altitude, a 20,000-ft/sec velocity, a one- inch nose radius, and a 
surface temperature of 2700°R. No mass injection and a noncatalytic 
wall are assumed. With the exception of the species profile, essentially 
no effect of diffusion model was observed. Even the effect on species 
profiles was small as shown in figures 27 and 28. 

Figure 27(a) and (b) show comparisons of binary and multicomponent 
species profiles for equilibrium air at an C of 0.0 and 1.0 respective- 
ly. More molecular nitrogen and less molecular oxygen are present at 
the wall for multicomponent diffusion. 

The elemental composition is not a constant for multicomponent 
diffusion as is the case for binary diffusion. However, the departure 
from elemental invariance is small for the equilibrium air calculation. 
The elemental composition at the wall was = 0.278 and = 0.722 



73 


at 5 = 0.0, and = 0.279 and = 0.721 at ? = 1.0, as opposed to 

the free stream values of C. = 0.24 and C„ = 0.76. 

0 N 

Figures 28(a) and (b) show comparisons of species profiles for 
multicomponent and binary diffusion models at ^ values of 0 and 1.0 
respectively. The flow chemistry is nonequilibrium. Once again the 
diffusion model has a small effect on species profiles. The elemental 
wall concentrations were = 0.221 and = 0.778 at ^ = 0 , and 
= 0.222 and = 0.779 at C = 1.0. 

The other flow parameters and surface transport values were 
essentially the same for both diffusion models, and therefore, compari- 
sons of results are not presented. However, table XIII, which presents 
dimensional heating rates for various body stations, illustrates this 
point. It can be seen from the table that the heat transfer rates for 
the flow conditions considered were about the same for both diffusion 
models. 


Mass Injection 

In the following, the results of computations with mass injection 
are considered. In all cases the mass injection rate distribution is 
given by equation (2.39). Figure 29 shows the injection distribution 
for different values of the nondimensional stagnation injection para- 
meter m^. It is seen that the injection rate approaches zero at a 
5 of about 1.5. 

Injecting air into air . -Figures 30 through 37 and table XIV show 
how injecting equilibrium air intb nonequilibrium air influences flow 
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parameters and surface transport. Multicomponent diffusion and a non- 
catalytic wall are assumed. Stagnation mass injection rates as large 
as 0.4 are considered. 

Figures 30(a) through (c) show comparisons of species profiles for 
no injection and an injection rate of 0.2. These comparisons are made 
at a ^ of 0., 1.0, and 3.0. At the stagnation station, the shock- 
layer chemical composition is altered substantially due to mass injec- 
tion (fig. 30(a)). At the wall, the mass of dissociated flow is 6 
percent for a 0.2 injection rate and 24 percent for no injection. 
Downstream the effect of mass . injection on. chemical composition decreases 
(fig. 30(b)) because the injection rate is smaller (eq. (2.39)). At 
3.0, which is beyond the body station where mass injection is zero, the 
chemical composition of' the two shock layers (fig. 30(c)) is about the 
same. 

Figures 31 and 32 show the stagnation tangential velocity and 
temperature profiles, respectively. These results qualitatively are as 
one would expect. That is, the shock- layer thickness increases and the 
tangential velocity and temperature gradients decrease with increasing 
mass injection rate. 

Figure 33 shows the effect of mass Injection on shock stand-off 
distance. For stagnation injection rates of 0.2 or greater, the 
injectant increases the shock stand-off distance as far downstream as 
the calculations are made (5 = 3.0). 

Figure 34 shows the effect of inejcting equilibrium air into non- 
equilibrium air on nondimens ional heat transfer distribution. As the 
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injection rate increases, the heat transfer rate near the stagnation 
point approaches zero for an of O.A. Downstream, however, the 

effect of injection rate on heat transfer rate decreases rapidly. These 

0 

results also show that the smaller injection rates are more effective 
in reducing the surface heating. Figure 35 shows that the same conclu- 
sions apply to Stanton number distributions. 

Figure 36 shows the effect of mass injection on skin friction 
coefficient distributions. In the mass injection region, large reduc- 
tions in skin friction occur. Downstream, however, the skin friction 
values with injection are approaching the no injection values. For 
example, for a stagnation mass injection rate of 0.4, the skin friction 
at ^ = 0.4 is 90 percent less than the no injection value, but at 
5 = 3.0 the reduction is only 9 percent. 

The effect of mass injection on wall pressure distribution is shown 
in figure 37. An almost negligible effect was obtained. The wall 
pressure is less with mass injection but the pressure differences 
never exceed 4 percent for the mass injection rates considered. This 
result is due primarily to the way that the injected mass alters the 
shock shape and, consequently, the shock pressure. This is evident 
in figure 33, where the shock slope decreases with increasing injection 
rate. Hence, the local shock angle, a, (fig. 1) is less and so is the 
shock pressure. Obviously, the way mass injection alters the shock 
shape and shock pressure will depend on the mass injection distribution 


which is assumed. 
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For some hypersonic flow conditions, in particular hypervelocity 
planetary entries where radiant heating is very large, the mass injec- 
tion rates can be large enough so that the convective heating is reduced 
to zero. Calculations were made to demonstrate the present capability 
to predict the effect of mass injection on stagnation heat transfer and 
shock shape at such injection rates. Figures 38 through 40 show such an 
effect for the injection of equilibrium air into reaching equilibrium air. 
Figure 38 shows how the stagnation wall heat transfer rate decreases with 
increasing mass injection — becoming zero at a mass injection rate of 
about 0.4. Figure 39 shows the influence of injection on stagnation 
shock stand-off distance. For the range of injection rates shown, the 
shock stand-off distance increases linearly with increasing injection 
rate. Figure 40 shows the effect of even larger injection rates on shock 
stand-off distance. (The results for injection rates greater than 0.4 
did not consider the downstream influence on the stagnation solution.) 

For injection rates greater than 0.4, the shock stand-off distance 
continues to increase with increasing injection rate but the slope of 
the curve decreases with increasing injection rate. 

Injecting water into air .-Injecting water either as a liquid or as 
a gas into a reacting flow is an example of transpiration cooling. A 
Qu^^ 0 ^t study (ref. 37 ) shows water to be one of the most efficient 
transpiration system fluids for limiting surface temperatures. This 
result was contingent on the transpiration cooling design being such 
that the water vaporizes at the wall-flow field interface. 
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Since water is an effective transpiration coolant, calculations 
were made to determine the effect of injecting water into reacting equi- 
librium air. The calculations were made for a 200,000 ft altitude, a 
20,000 ft /sec velocity, a 45° hyperboloid with a 12-inch nose radius, 
and a wall temperature of 2700°R. The large nose radius (12 inches 
vs 1 inch) was used so that the equilibrium flow approximation would be 
more realistic. That is, previous results have shown that the equili- 
brium flow approximation is not valid for the above free stream conditions 
and a one-inch nose radius calculation. 

The Injected water (gas phase) was in equilibrium with the wall 
temperature and pressure and consisted of a total of six possible 
chemical species. These species were H, 0, ^ 2 > *^ 2 * H 2 O. 

Figures 41 (a) through (c) show the equilibrium composition of water for 
the aforementioned six species over a temperature range of 1800° to 5400° 
R and for pressures of 0.1 and 0.01 atmosphere, respectively. Since 
the calculations were for a wall temperature of 2700°R and wall pressures 
of 0.1 atmosphere or less, the injectant was essentially all H^O. 

(See figure 41(a)). 

The equilibrium flow is assumed to consist of nine species, the six 
previously mentioned plus N, N 2 , and NO. Since the larger number of 
chemical species increases the computational requirements, downstream 
solutions are presented only as far as ? = 1.6, and the solutions are 
iterated in 5 only once. Also the diffusion is binary where the 
diffusion coefficient is taken as the coefficient for molecular nitrogen 
diffusing into atomic oxygen. 
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Figure 42 shows how the stagnation elemental profiles are influenced 
by H 2 O (gas) injection rates. Note that for an injection rate of 0.1, 
the elemental hydrogen mass fraction at the wall is 0.11, the same as 

that for H^O. For the 0.05 injection rate, the hydrogen elemental mass 

fraction at the wall is 0.095. Figure 42 also shows the distance that 
the elemental hydrogen is transported away from the wall. The elemental 

hydrogen mass fraction is 0.001 at n = 0.24(n = 0.014) and at 

n = 0.31 (n = 0.019) for 0.05 and 0.1 injection rates, respectively. 

Species profiles for H 2 O injection rate of 0.1 are shown in figures 
43(a) through (c) for C values of 0, 0.4, and 1.0. Figures 44(a) 
through (c) show the species profiles for a 0.05 injection rate at the 
same S values. For the range of C considered, the wall mass 
fraction of H 2 O is approximately 1.0 for the 0.1 injection rate and 
approximately 0.85 for the 0.05 injection rate. (A detailed listing of 
the stagnation wall species mass fractions is included in table XV.) 
Dissociation occurs as the H 2 O (gas) is transported into the higher 
temperature flow. The dissociation of H 2 O is strongly endothermic and 
substantially reduces the flow field temperature. 

Figure 45 shows the stagnation temperature profiles for H 2 O injection 
rates of 0 .0 , 0 .05 , and 0 . 1. At the wall, the temperature gradients 
decrease with increasing injection rate. As a result of the H 2 O 
injection, the temperature profiles experience a marked reduction in 
temperature for about one-third of the shock-layer thickness adjacent 
to the wall. This temperature reduction occurs because energy is 
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absorbed as the sensible enthalpy of the H 2 O is increased and as the 
H 2 O undergoes endothermal reactions. 

Dissociation of H 2 O is highly endothermic. For example, the dissocia- 
tion of one gram of H 2 O at 4500°R (2500°K) , forming H and 0, requires 
12,7 kcal/gram. In contrast, the dissociation of one gram of O 2 at 
the same temperature requires only 3. 8 kcal/gram. 

Figure 46 shows the effect of H 2 O injection on stagnation tantential 
velocity profiles. Note that the velocity profiles have a distinct 
boundary-layer region. The effects of H 2 O injection on heat transfer, 
Stanton number and skin friction coefficient are compared with the 
results for air and ablation injection in a subsequent section. 

Injecting ablation products into air . -When an ablator is used to 
protect the surface of a vehicle from a high- temperature environment, 
gaseous species are injected into the flow field as the ablator is 
thermally degraded. The ablative products injected into the flow field 
have the same qualitative effect as that obtained with transpiration 
cooling. For an ablating surface, however, the number of gaseous 
species injected into the flow field is generally much larger than that 
for a transpiration cooled surface. Consequently, analyses for reacting 
flows adjacent to ablating surfaces must account for a large number of 
reacting species. 

To determine the effect of injecting ablation species on flow 
parameters and surface transport, calculations were made for a low- 
density phenolic-nylon ablator. Gaseous ablation products in chemical 
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equilibrium at the wall temperature and pressure were injected into 
equilibrium air. Results are presented for Injection rates of 0.0, 

0.05, and 0.1. The calculations were made for the same free stream and 
wall temperature conditions as that for the H^O injection calculations. 

The elemental composition used for the low— density phenolic nylon 
ablator is that given by reference 38 and is = 0.699, = 0.163, 

= 0.0812 and = 0.0568. For equilibrium calculations, the elemental 
composition completely specifies the material. A total of nineteen 
possible species are assumed to be injected into the equilibrium flow 
field. The species are those listed in table II with the exception of 
OH and H 2 O. The equilibrium species mass fractions for the aforementioned 
elemental composition and assumed species are shown as a function of 


temperature in figures 47(a) and (b) for pressures of 0.1 and 0.01 
atm, respectively. 

The possible flow field species are assumed to be the same as those 
ax ;0 injected at the wall. For this system of species, the computa- 
tional time required to obtain solutions become quite large as indicated 
in table VI. Consequently, downstream solutions are presented only 
as far as C = 0 . 6 , and the solutions are iterated in K only once. 

The diffusion is binary where the diffusion coefficient is taken as the 
coefficient for molecular nitrogen diffusing into atomic oxygen. 

Since the wall temperature is 2700°R (1500°K) , the injected ablation 
species are predominantly ^-^2 shown in figures 47(a) and 

(b) . Also, the species at the wall are essentially C 2 H 2 and CO. A 
detailed listing of the stagnation species mass fractions and other wall 
properties are included in table XV. 
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Figure 48 presents the stagnation elemental profiles for ablation 
injection rates of 0.05 and 0.1. For these injection rates, the carbon 
and hydrogen elements are confined to about one third or less of the 
shock-layer thickness adjacent to the wall. 

For the 0.1 injection rate, figures 49(a) through (b) show the species 
profiles at ^ values of 0.0 and 0.4. The species profiles for the 
same ^ values and a mass injection rate of 0.05 are shown in figures 
50(a) and (b) . These results show that a large number of chemical 
reactions occur near the wall. 

The chemical reactions that remove C 2 H 2 are rather important because 
the chemical process is highly endothermic and C 2 H 2 is the major species 
at the wall. At the stagnation point, the mass fraction of C 2 H 2 is 
0.574 for a 0.1 injection rate and 0.446 for a 0.05 injection rate. 
However, the mass fraction of C 2 H 2 decreases rapidly with increasing 
n. At an n of 0.15, the mass fraction is 0.015 for the 0.1 injection 
rate and zero for the 0.05 injection rate. 

The second most abundant species at the wall is CO. For the tempera- 
ture interval in which the carbon element is available (fig. 48), CO 
dissociation does not occur (see figures 47(a) and (b)). In fact, addi- 
tional CO is formed as some of the hydrocarbons are chemically degraded 
at the higher temperatures . 

The effect of injecting ablation products on stagnation temperatures 
is shown in figure 51. Results for 0.0, 0.05, and 0.1 injection rates 
are shown. Near the wall, the temperatures are substantially reduced by 
injection of ablation products. For example, 35 percent of the shock 
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layer adjacent to the wall experiences a temperature reduction for a 0.1 
injection rate and 30 percent of the shock, layer experiences a temperature 
reduction for a 0.05 injection rate. 

The effects of injecting ablation species on stagnation velocity 
profiles are shown in figure 52. The effect of ablation injection on 
heat transfer, Stanton number, and skin friction coefficient is compared 
in the following section. 

Comparison of results for mass injection of air, H^O (gas), and 
abalation species . -In this section, comparisons of the effects of 
injecting air, H 2 O, and ablation species on heat transfer rate, Stanton 
number, and skin friction distributions are presented. The equilibrium 
air calculations were made by using the same free stream and wall 
temperature conditions stated in the two previous sections. For the 
equilibrium air results, multicomponent diffusion was used, and the 
solutions were iterated twice in 5. (See table XVI for a listing of 
wall pressure ratio, shock stand-off distance, and the aforementioned 
quantities at K values of 0, 0.2, 0.4 and 0.6.) 

Figure 53 shows the heat transfer rate distributions for injection 
rates of 0.0, 0.05, and 0.1. These results show that the mass injec- 
tion in the stagnation region of any of the three injectants significant- 
ly reduces the heat transfer rates. Downstream of the mass injection 
region (see fig, 53), the H 2 O injectant provides a larger reduction in 
heating than the air injectant (for ablation Injection, the solutions 
presented are only within the injection region). The heat transfer 
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rates with air injection approach the no injection values a short 
distance down stream of the injection region. This agrees with the 
previously discussed effect of air injection for a one— inch nose radius 
body. For the 0.05 injection rate, H^O effects the largest reduction 
in heating followed by ablation products and finally air. For the 0.1 
injection rate, the H 2 O and ablation injectants effect about the same 
reduction in heating and both Injectants give more reduction than air 
injection . 

Figure 54 shows the effect of each of the three injectants on 
Stanton number distribution. Most of the conclusions concerning the 
effect of mass injection on heat transfer distribution are also appli- 
cable to Stanton number distribution. However, the relative effective- 
ness of the three injectants in reducing Stanton number is not the same 
for the 0.05 injection rate as that for reducing the heat transfer. 

This effect is due to large variations in wall enthalpy (see table XV) 
for the three injectants. 

Figure 55 shows the effect of mass injection on skin friction 
coefficient and how this effect differs for air, H 2 O and ablation 
injection. Large skin friction reductions occur in the mass injection 
region. Downstream of the injection region, the skin friction values 
approach the no injection values. In the mass injection region, the 
three injectants effect about the same reduction in skin friction. 
Downstream of the injection region, the H 2 O injectant provides a larger 
reduction in skin friction than does the air injectant. 
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VII. CONCLUDING REMARKS 

Equations describing reacting and nonreacting viscous shock layers 
about an axisymmetric body at zero angle of attack were presented. 

These equations account for mass injection and multicomponent diffusion. 
An implicit finite-difference technique for solving these equations was 
discussed. Calculated results were described which show the effect of 
chemistry model, surface catalyticity , diffusion model, and mass injec- 
tion on flow parameters and surface transport. Some of the solutions 
provide information previously unavailable for reacting air. For 
example, the downstream viscous shock-layer solutions that treated the 
diffusion as mutlicomponent, with or without injection of air, were 
previously unavailable. These results were obtained for both non- 
equilibrium and equilibrium flow. 

Results of the study lead to the following conclusions: 

1. The chemistry model substantially influences flow parameters 
and surface transport. Results show that the frozen and equilibrium 
chemistry models overpredict wall heat transfer rates for the flow 
conditions considered in this study. For example, the stagnation heating 
rate obtained with an equilibrium calculation was 54 percent greater 
than comparable results with a nonequilibrium and noncatalytic wall 
calculation. 

2. Wall catalyticity significantly influences wall heat transfer 
rate and species composition near the wall but has little influence 

on the other flow parameters considered. For an equilibrium catalytic 
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wall, the wall heat transfer rate for nonequilibrium flow was shown to 
be about the same as that for equilibrium flow. 

3. For reacting air, results were compared for a multicomponent 
and for a binary diffusion model. For the binary approximation, the 
diffusion coefficient was the diffusion coefficient for atomic oxygen 
diffusing into molecular nitrogen. Results of the comparisons show 
negligible differences in most flow parameters and surface transport 

for the two diffusion models. This result applies to both nonequilibrium 
and equilibrium air. For multicomponent diffusion a small amount of 
elemental diffusional separation was observed, and consequently, a 
small difference in species profiles occurred. 

4. The effect that mass injection has on heat transfer rates de- 
creases as the injection rate increases. Also, the effect that mass 
injected in the stagnation region has on the flow field decreases very 
rapidly downstream. 

Results of stagnation solutions were presented for mass injection 
rates substantially in excess of those required to reduce the convective 
heating rateto zero. For the injectants considered, water was most 
effective in reducing wall heat transfer rates. 

5. The shape of the equilibrium temperature profiles differs from 
the results of a previous stagnation viscous shock-layer solution; 
however, the shape of the present temperature profiles can be sub- 
stantiated by considering the chemical reactions that occur within the 


shock layer. 
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6. Two areas for improving the present analyses are noteworthy. 
First, it is desirable to reduce the time required for the equilibrium 
flow calculations. This can be accomplished by developing Improved 
equilibrium chemistry computer codes and by developing more efficient 
means for solving the energy equation. Secondly, the range of shock 
Reynolds numbers for which the nonequilibrium analysis can be applied 
needs to be extended. All the nonequilibrium results reported were 
for a shock Reynolds number of 191. If the shock Reynolds number is 
increased beyond about 500, the analysis will not converge. 

7, These analyses should be appropriate for calculating radiating 
viscous flows with mass injection. 
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TABLE IV. -REACTION RATE COEFFICIENTS 















TABLE V.-FREE STREAM CONDITIONS 
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TABLE VI. -SUMMARY OF COMPUTING TIME AND ITERATION 
REQUIREMENTS FOR VARIOUS CASES^ 
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TABLE VII. FROZEN AIR STAGNATION SHOCK AND WALL VALUES 
[alt = 200,000 ft; U* = 20,000 ft/sec; a* = 1.0 in.; 

Re = 2.597X103; Re = I.OIAXIO^] 
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TABLE VIII. -EQUILIBRIUM AIR STAGNATION SHOCK A^JD WALL VALUES 
[alt = 200,000 ft; U* = 20,000- ft/sec; a* = 1.0 in.; Re^ = 
2.597X103; Re = 3.376X10^] 
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TABLE IX. -EQUILIBRIUM AIR STAGNATION SHOCK AND WALL VALUES^ 

[alt = 280,000 ft; U* = 26,000 ft/sec; a* = 12 in.; Re^ = 

1.482X10 3; Re = 1.498X102 
s 


PROPERTIES 

SHOCK 

WALL 

Temperature, °R 

1.032X10^ 

2.429X10^ 

Pressure, atm , . . . 

4.396X10"^ 

4.528X10"® 

3 

Density, Ibm/ f t 

9.274X10"^ 

7.411X10"® 

2 


. -6 

Viscosity, Ibf-sec/ft 

2.513X10 

1.032X10 

Velocity, ft/sec 

-1.306X10^ 

0.0 

Enthalpy, Btu/lbm 

1.361X10^ 

5.061X10® 

Molecular weight, g/g-mole. . . . 

1.580X10^ 

2.888X10^ 

Lewis number 

1.400X10® 

1.400X10® 

Prandtl number 

6.734X10"^ 

6.957X10"^ 

Mass fractions: 

0 

2.398X10"^ 

5.331X10"® 

0. 

2.098X10"® 

2.396X10"^ 

1 

N 

5.927X10"^ 

5.293X10"®® 

N„ 

1.671X10"^ 

7.597X10"® 

1 

NO 

3.776X10"^ 

6.008X10”'^ 


= 3.801X10 q* = 6.927X10^ Btu/ft^-sec, 
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TABLE X. -EQUILIBRIUM AIR STAGNATION SHOCK AND WALL VALUES® 
[alt = 150,000 ft; U* = 22,520 ft/sec; a* = 0.5 in.; Re^ = 
9.236X10^; Re = 1.134X10-^ 


PROPERTIES 

SHOCK 

L 

WALL 

Temperature, °R 

1.221X10^ 

-1 

2.657X10^ 

-1 

Pressure, atm 

7.729X10 

8.038X10 

3 

Density, Ibm/ft . • * . . • • * • 

- 

1.619X10 

1.197X10'^ 

2 

Viscosity, Ibf-sec/ft 

2.861X10"^ 

1.094X10"^ 

Velocity, ft/sec 

-1.547X10^ 

o 

o 

Enthalpy, Btu/lbm , 

1.020X10 

5.411X10 

Molecular weight, g/g-mole. . . . 

1.863X10 

2.888X10 

Lewis number 

1. 400X10*^ 
~1 

1.400X10° 

Prandtl number 

Mass fractions: 

6.627X10 

6.958X10 ^ 

0 

2.371X10”^ 

-4 

3.638X10"^ 

0^ 

1.108X10 

2.397X10 ^ 

2 

N 

3.260X10'^ 

-1 

3.983X10"^^ 

-1 



4.310X10 

7.593X10 

2 

-3 

-4 

NO . 

J 

5.712X10 

9.205X10 


= 5.601X10“^; q* = 2.375X10^ Btu/ft^-sec. 




TABLE XI. -NONEQUILIBRIUM AIR STAGNATION SHOCK AND WALL VALUES 
[alt = 200,000 ft; U* = 20,000 ft/sec; a* = 1.0 in.; Re^ = 
2.597X103; Re = 1.914X10^ 
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TABLE XII. -EQUILIBRIUM AIR STAGNATION SHOCK AND WALL VALUE 
[alt = 100,000 ft; U* = 20,000 ft/sec; a* = 1.0 in.; Re^ = 
1.791X105; Re = 1.876X10'^] 


Properties 

Shock 

Wall 

Temperature, °R 

1.255X10^ 

2.700X10° 

6.031X10° 

Pressure, atm 

5.771X10° 

3 

Density, Ibm/ft . . . . 

1.312X10'^ 

8.913X10 

2 

Viscosity, Ibf-sec/ft 

2.948X10"^ 

I.IIOXIO"^ 

Velocity, ft/sec 

-1.627X10° 

0.0 

Enthalpy, Btu/lbm 

8.019X10° 

5.79 7X10° 

Molecular, weight, g/g-mole . . . . 

2.897X10^ 

2.076X10^ 

Lewis number^ 

- 

- 

Prandtl number 

6.598X10"^ 

6.964X10"^ 

Mass fractions: 

0 

2. 305X10"^ 

7.256X10"® 



6.336X10”'^ 

2.640X10"^ 

2 


-14 

N . . . 

1.778X10 

1.044X10 



5.746X10"^ 

7.346X10"°- 

2 

-2 

-3 

NO 

1.658X10 

1.384X10 


= 6.708X10 q* = 2.735X10^ Btu/ft^-sec. 


b 


Multicomponent diffusion. 




TABLE XIII. -SUMMARY OF DIMENSIONAL HEAT TRANSFER RATES AT DIFFERENT BODY STATIONS 
[alt = 200,000 ft; U* = 20,000 ft/sec; a* = 1.0 in.; Re = 2.597X10^] 
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TABLE XIV. -NONEQUILIBRIUM STAGNATION WALL CONDITIONS FOR A NONCATALYTIC WALL AND MASS INJECTION OF 
EQUILIBRIUM AIR [alt = 200,000 ft; U* = 20,000 ft/sec; a* = 1.0 in.; Re^ = 2.597X103; Re = 

1.914X10'^; multicomponent diffusion] 
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See TABLE XI for shock properties. 












TABLE XV. -ST AGNATION EQUILIBRIUM WALL CONDITIONS FOR MASS 
INJECTION OF WATER AND ABLATION SPECIES [alt = 200.000 
ft: U* = 20,000 ft/sec; a* = 12 in.; Re = 3.116X10^; Re 

*00’ ’o ^ CO S 

2.297X10^; binary diffusion] 
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TABLE XVI.-NONDIMENSIONAL HEAT TRAfISFER-SKIN FRICTION COEFFICIENT, STANTON 
NUMBER, WALL PRESSURE RATIO, AND SHOCK STAND-OFF DISTANCE VALUES FOR NO IN- 
JECTION AND INJECTION OF AIR^, WATER^ , AND ABLATION PRODUCTS^ AT DIFFERENT 

BODY STATIONS 
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Multicomponent diffusion 
Binary diffusion 



TABLE XVI. -NONDIMENS ION AL HEAT TRANSFER SKIN FRICTION COEFFICIENT, STANTON 
NUMBER, WALL PRESSURE RATIO, AND SHOCK STAND-OFF DISTANCE VALUES FOR NO IN- 
JECTION AND INJECTION OF AIR, WATER, AND ABLATION PRODUCTS AT DIFFERENT BODY 

STATIONS - Concluded. 
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PvCi Ji 



^ Body surface 

(pvCi)_ 

(a) Mass balance at gas- solid interface. 



(b) Energy balance at gas- solid interface. 

i 

Figure 2.- Surface boundary conditions for mass and energy transfer. 






















Guess a value for Pg 


Solve Eq. (2.49d) for (hg)i 



* 

Solve Eq. (2.49c) for Pg 


* 

Solve Eq. (2.49e) for Tg 




(b) Shock frozen at free stream chemistry. 


Figure 3.- Concluded. 
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(a) Equilibrium chemistry. 


Figure 4.- Flow chart for solution sequence of viscous shock-layer 

equations. 
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(b) Skin friction variation with distance along body surface. 
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(b) Skin friction variation with distance along body surface. 





(c) Heat transfer variation with distance along body surface. 
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Temperature, °R 


(a) Stagnation temperature profiles for alt * 280,000 ft, U*oo = 26,000 
ft/sec, and a* * 12 in. 


Figure 7.- Comparison of the present calculations with the stagnation 
results of Edelman and Hoffman for equilibrium air. 
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(c) Stagnation temperature profiles for alt ■ 150,000 ft, U*«o « 22,520 
ft/sec, and a* = 0.5 in. 


Figure 7.- Continued, 
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Temperature, °R 


(a) Stagnation temperature profiles. 

Figure 8,- Comparison of the present calculations with the stagnation 
results of Blottner for nonequilibrium air. 
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Figure 8,- Continued. 
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(c) Stagnation species profiles. 


Figure 8,- Concluded. 
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(c) Frozen chemistry. 
Figure 9. - Concluded. 
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(a) s = 0.0. 


Figure 11.- Species profiles for nonequillbrium air with multi-component 

diffusion and no mass injection. 
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(a) Temperature profile. 


Figure 14.- Stagnation equilibrium air results for an alt = 100,000 ft 

U* - 20,000 ft/sec, and a* = 1,0 in. 
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(a) s = 0.0. 


Figure 15.- Effect of chemistry model on tangential velocity profiles with 
multicomponent diffusion and no mass injection. 
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(b) s = 2.0. 


Figure 15.- Concluded 
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Figure 17.- Effect of chemistry model on shock stand-off distance with 
multicomponent diffusion and no mass injection. 
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Figure 18. Effect of chemistry model on nondimensional heat transfer with 
multicomponent diffusion and no mass injection. 
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Figure 19.- Effect of chemistry model on Stanton number with multicomponent 

diffusion and no mass injection. 
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Figure 20.- Effect of chemistry model on skin friction coefficient with 
mult X component dxffusxon end no mess injection. 
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Figure 21.- Effect of chemistry model on wall pressure distribution with 
multicomponent diffusion and no mass injection. 
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Noncatalytic wall 
Equilibrium catalytic wall 

alt = 200,000 ft 
U* = 20,000 ft/sec 
T* = 2700° R 
a* = 1.0 in. 

Multicomponent diffusion 
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- Equilibrium catalytic wall 
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Multicomponent diffusion 






Figure 23.- Effect of wall catalytic! ty on shock stand-off distance 
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Figure 24.- Effect of wall catalyticity on nondimens ional heat transfer. 
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Figure 27.- Effect of diffusion model on species profiles for equilibrium 






Mass fractions, 


O 



n/ng 

(a) s = 0. 


Figure 28.- Effect of diffusion model on species profiles for nonequilibrium 
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(a) s = 0.0. 
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Figure 32.- Stagnation temperature profiles for different injection rates 
of equilibrium air into reacting nonequilibrium air. 
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Figure 33.- Comparison of shock stand-off distances for different injection 
rates of equilibrium air into reacting nonequilibrium air. 
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Figure 34.- Comparison of nondimensional heat transfer distributions for 
different injection rates of equilibrium air into 
reacting nonequilibrium air. 
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Figure 35.- Comparison of Stanton number distributions for different in' 
jection rates of equilibrium air into reacting non- 
equilibrium air. 
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Figure 36.- Comparison of skin friction coefficient distributions for 
different injection rates of equilibrium air into 
reacting non equilibrium air. 





Figure 37.- Comparison of wall pressure ratio distributions for different 
injection rates of equilibrium air into reacting 
nonequilibrium air. 
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Figure 38.- Stagnation heat transfer as influenced by injecting equili 

brium air into reacting equilibrium air. 




Figure 40.- Stagnation shock stand-off distance as influenced by large mass 
injection rates of equilibrium air into reacting equilibrium 
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Temperature, °K 
(b) Pressure = 0.01 atm. 

Figure 41.- Concluded. 
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Figure 44.- Continued 
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Figure 44.- Concluded. 
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Figure 45.- Stagnation temperature profiles for different injection rates 

of H„0 into equilibrium air. 
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Figure 46.- Stagnation velocity profiles for different injection rates of 

H_0 into equilibrium air. 
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Figure 48.- Stagnation elemental profiles for Injecting ablation species 

into equilibrium air. 
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Figure 51.- Stagnation temperature profiles for different injection rates 

of ablation products. 
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Figure 52.- Stagnation velocity profiles for different injection rates of 

ablation products. 
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Figure 53.- Comparison of nondlmenslonal heat transfer distributions for 
different Injection rates of air, H2O, and ablation 
species Into equilibrium air. 
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Figure 54.- Comparison of Stanton number distributions for different injec- 
tion rates of air, H 2 O, and ablation species into 

equilibrium air. 
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Figure 55.- Comparison of skin friction coefficient distributions for 
different injection rates of air, H^O, and ablation 
species into equilibrium air. 
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